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NOMENCLATURE 


A  Cross-sectional  area 

d  Distance 

C  Specific  heat 

Evaporation 
m  Mass  flowrate 
M  Molecular  weight 

P  Pressure 

Q  Heat  transfer 

R  Universal  gas  constant 

r  Radius  of  curvature 

t  Time 

T  Temperature 
V  Velocity 

W  Wetted  perimeter 

X  Length  or  distance 

y  Distance 

Greek  Alphabet: 

3  Dimensionless  parameter  (Fig.  2) 
e  Wetting  angle 
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X  Latent  heat 
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Subscripts: 
a  Adiabatic 

b  Boundary 

c  Capillary,  curvature 

d  Diagonal 

e  Evaporator 
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1  Liquid 

m  Momentum,  mean 

o  Outer 

P  Pipe 
s  Solid 
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1.0  INTRODUCTION 


The  need  to  develop  microelectronic  devices  capable  of  operating  at 
increased  performance  levels  with  high  reliability,  requires  better  methods  of 
thermal  control.^  With  the  recent  trend  towards  increased  miniaturization  and 
component  density,  the  thermal  management  problem  increases  dramatically.  This 
is  particularly  true  in  the  microelectronic  components  proposed  for  use  in 
spacecraft  applications  where  the  physical  space  is  at  a  premium,  and  the  size 
of  the  heat  removal  system  is  a  major  contributor  in  the  determination  of 

overall  spacecraft  size. 

2 

In  1984,  Cotter  described  very  small,  miniature  heat  pipes  for  cooling 
microelectronic  devices.  He  defined  a  micro  heat  pipe  as  "one  so  small  that  the 
mean  curvature  of  the  vapor-liquid  interface  is  necessarily  comparable  in 
magnitude  to  the  reciprocal  of  the  hydraulic  radius  of  the  total  flow  channel." 
In  practical  terms,  a  micro  heat  pipe  is  a  wickless,  noncircular  channel  with  an 
approximate  diameter  of  100-1000  ym  and  a  length  of  about  10-20  mm.  These  micro 
heat  pipes  do  not  have  a  wick  but  use  the  mean  curvature  of  the  channel  to  pump 
the  liquid  from  the  condenser  to  the  evaporator.  In  addition  to  the 
difficulties  associated  with  the  manufacturing,  cleaning,  and  filling  of  these 
devices,  they  are  very  sensitive  to  the  shape,  size,  and  amount  of  working  fluid 
present.  Since  the  liquid-vapor  interface  changes  continually  along  the  pipe, 
care  must  be  exercised  to  ensure  proper  wetted  conditions  without  flooding. 

Under  a  previous  contract  with  Wright-Pat terson  AFB,  steady-state 
experimental  and  analytical  investigations  of  a  trapezoidal  micro  heat  pipe  were 
conducted  to  identify  and  better  understand  the  phenomena  which  govern  the 
performance  limitations  and  operating  characteristics.  The  analytical  portion 
of  the  investigation  began  with  the  development  of  a  steady-state  model,  in 
which  the  effects  of  the  extremely  small  characteristic  dimensions  on  the 
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conventional  steady-state  heat  pipe  modeling  techniques  were  examined. 

In  addition  to  the  modeling  portion  of  the  investigation,  several  micro 

2 

heat  pipes  approximately  1  mm  in  cross-sectional  area  and  57  to  60  mm  in 

length,  were  evaluated  experimentally  to  determine  the  accuracy  of  the  model  and 

provide  verification  of  the  micro  heat  pipe  concept.  Tests  were  conducted  in  a 

vacuum  environment  to  eliminate  conduction  and  convection  losses.  The 

steady-state  experimental  results  obtained  were  compared  with  the  analytical  and 

numerical  models  and  verified  the  modeling  predictions  for  both  the  steady-state 

model  and  the  transient  model  at  steady-state.  The  results  indicated  that  the 

vetting  angle  was  the  single  most  Important  factor  contributing  to  the  transport 

capacity  and  thermal  behavior  of  these  types  of  heat  pipes. 

Micro  heat  pipes  have  a  wide  range  of  applications.  Of  particular  interest 

here  is  the  thermal  control  of  discrete  semiconductor  devices.  Other 

applications  include  the  thermal  control  of  infrared  detectors  and  the  cooling 

of  the  leading  edges  of  hypersonic  aircraft.  Smaller  micro  heat  pipes,  (less 
2 

than  0.05  mm  in  cross-sectional  area),  could  be  incorporated  into  semiconductor 
devices  as  an  integral  part  of  the  chip.  This  would  provide  a  mechanism  for 
eliminating  the  localized  hot  spots  associated  with  semiconductor  devices  and 
significantly  improve  the  chip  reliability.  In  addition,  these  micro  heat  pipes 
could  provide  a  method  for  collecting  and  dissipating  the  high  thermal  loads 
associated  with  the  three-dimensional  chip  architectures  being  developed  for  use 

3 

in  fifth  generation  computers. 

1.1  Objectives  of  the  Investigation 

Although  the  models  developed  previously  were  shown  to  accurately  predict 
the  steady-state  performance  limitations  and  operational  characteristics,^  no 
experimental  data  on  the  transient  operation  and  characteristics  were  obtained. 
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In  addition,  several  factors  which  significantly  affect  both  the  steady-state 
and  transient  behavior  required  continued  study.  To  further  understand  the 
behavior  of  these  types  of  devices,  a  transient  model  was  developed  to  predict 
the  behavior  during  startup  or  variations  in  the  evaporator  thermal  load. 

This  follow-on  investigation  included  several  major  tasks.  The  first  of 
these  involved  the  expansion  of  the  existing  analytical  model  to  include  the 
capability  of  accurately  predicting  the  transient  performance  of  micro  heat 
pipes.  Second,  was  a  determination  of  the  effect  of  variations  in  the 
cross-sectional  shape  of  the  device  and  the  optimization  of  this  shape.  Third, 
was  the  experimental  investigation  of  a  micro  heat  pipe  which  utilized  this 
optimum  shape  to  remove  heat  from  a  ceramic  chip  carrier.  And  finally,  was  a 
comparlslon  of  the  analytical  predictions  and  the  experimental  results. 

The  following  report  is  divided  into  two  major  sections:  the  first  deals 
with  a  description  of  the  development  and  expansion  of  the  transient  numerical 
model,  while  the  second  describes  the  experimental  program  and  compares  the 
modeling  predictions  and  the  experimental  results. 


2.0  DESCRIPTION  OF  THE  NUMERICAL  MODEL 

4 

Babin  et  al.  described  a  micro  heat  pipe  as  one  which  satisfies  the 
condition 

^>1 

where  r^  is  the  capillary  radius  of  the  heat  pipe  and  r^^  is  the  hydraulic  radius 
of  the  flow  channel.  In  these,  as  in  conventional  heat  pipes,  the  high  thermal 
conductivity  is  the  result  of  the  evaporation  and  condensation  process  occurring 
within  the  heat  pipe.^'^  For  this  reason,  determination  of  the  evaporation  and 
condensation  rate  of  the  liquid  plays  a  key  role  in  evaluating  the  thermal 


3 


characteristics  and  heat  transport  limitations  of  these  heat  pipes.  In  the 

present  model,  an  expression  for  the  molecular  flux  of  evaporation,  J,  presented 

7  8 

by  Collier  and  later  used  by  Colwell  and  Chang  was  employed.  This 

relationship  is 


j  = 


M 


I  2ltiRT  J 


-  fv) 


(2) 


where  M  is  the  molecular  weight  of  the  fluid,  R  is  the  universal  gas  constant,  T 

is  the  absolute  temperature,  is  the  saturation  pressure,  and  is  vapor 

pressure.  When  P  is  greater  than  P  ,  j  is  positive  and  the  liquid  evaporates, 
s  ▼ 

When  P^  is  less  than  P^,  j  is  negative  and  the  vapor  condenses.  In  the 
development  of  the  numerical  model  presented  here,  the  evaporation- condensation 
rate  was  assumed  to  be  proportional  to  the  liquid-vapor  interface  area  in  each 
section  of  the  heat  pipe,  i.e. 


dM 


ec 


=jW^^dx 
e  ^  ec 


(3) 


where  is  the  evaporation-condensation  surface  area  in  one  section,  W^^  was 

the  perimeter  of  the  liquid-vapor  interface,  and  Ax  is  the  section  length. 

For  a  given  latent  heat,  X,  the  rate  of  the  heat  released  or  absorbed  in  any 
section  can  be  determined  from: 


AQ  =  XAM 
ec  ec 


(A) 


In  order  for  a  heat  pipe  to  operate  properly,  the  capillary  pressure  must  be 
sufficient  to  overcome  the  difference  of  the  liquid  and  vapor  pressures,  or 


P  ,  =  P  -  P,  =  a/r 
vl  V  1  c 


(5) 


where,  P^  is  the  vapor  pressure,  a  is  the  liquid  surface  tension  and  r^  is  the 
radius  of  curvature  of  the  liquid-vapor  interface.  Since  both  P^j^  and  r^  are 
functions  of  the  axial  position,  the  derivative  of  the  above  equation  yields: 
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(6) 


dP,,  -a  dr 

vl  _ _ c 

dx  r  ^  dx 

c 

Expressing  this  over  a  finite  interval  yields 


where 


>  T  =  - TT  Ar 

\)1  2  c 

r 

c 


vl  fv  mv  fl  ml  g 

The  pressure  changes  are  those  due  to  the  liquid  and  vapor  friction  and  the 
momentum  changes  are  those  occurring  in  any  one  section  along  the  heat  pipe. 
Based  on  the  momentum  equation,  the  pressure  drops  due  to  friction  for 


laminar  flow  can  be  expressed  as: 


2yAxV 


where  u  is  the  viscosity  of  the  fluid,  rj^  is  the  hydraulic  radius,  and  V  is  the 
velocity  of  the  flow  which  is  equal  to  a  function  of  the  mass  flow  rate,  the 
density,  and  the  area. 


The  total  pressure  drop  due  to  the  momentum  change  can  be  expressed  as: 


Vfiin  +  mAV  +  AraAV 


The  energy  equation  used  in  this  model  can  be  written  as 


where 


PlAiCpiATi( 


Q  „  +  Q 
out  r 


p  A„Co  AT  dx 
s  s  Ps  s 


Cpi*  Cp^  and  Cp^  are  the  heat  capacities  of  the  liquid,  vapor  and  solid 
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boundaries  respectively  and  the  AT's  are  the  corresponding  temperature  changes. 


The  continuity  equation  can  be  written  as 


Pv 


dx  =  constant  =  fluid  charge 


(13) 


and  the  total  volume  inside  the  heat  pipe  can  be  defined  and  expressed  as: 

Vol  =  Volj^  +  Vol^  (14a) 

and 

A  =  (14b) 

respectively. 

As  is  the  case  for  larger  heat  pipes,  micro  heat  pipes  consist  of  three 
basic  sections,  an  evaporator  and  a  condenser,  separated  by  an  adiabatic 
section.  Each  of  these  sections  has  a  different  set  of  boundary  conditions,  and 
as  a  result,  have  been  treated  independently. 

2.1  Evaporator  Region 

The  single  boundary  condition  utilized  in  the  evaporator  section  is  the 
time  dependent  heat  flux.  For  a  specific  input  heat  flux,  the  saturation 
pressure  at  a  given  location  can  be  obtained  by  using  a  combination  of  the 
molecular  flux  expression  given  in  Eq.  (2)  and  the  energy  conservation  equation 
given  in  Eq.  (12).  Assuming  that  a  temperature  change  at  the  boundary  yields  a 
corresponding  temperature  change  of  the  liquid,  these  equations  take  the 
following  form: 


P 

s 


AQ 


ec 


(M/2nRT^)^'^^  AxW 


+  P 


AQi  = 


ec 

Cp  Pj^Aj^  AxAT^ 

At 


AQ, 


Cn  pA  Ax  AT 
Ps  s  s  s 

At 


(15) 

(16) 

(17) 


6 


Since  Tj^  is  a  function  of  P^,  Eqs.  15  through  18  are  coupled  and  can  be 
solved  using  an  iterative  method  with  relaxation  to  obtain  values  for  ^ec’ 
CQ^,  and  Because  the  difference  between  the  boundary  and  liquid 

temperatures,  is  proportional  to  the  input  heat  flux,  the  boundary 

temperature  can  be  obtained  by  adding  to 

2.2  Adiabatic  Region 

In  the  adiabatic  section,  the  latent  heat  absorbed  by  evaporation  and  given 
up  by  condensation  is  equal  to  the  sensible  heat  absorbed  or  rejected  by  the 
liquid  and/or  the  heat  pipe  case.  The  governing  equations  for  this  section 
can  therefore  be  expressed  as: 


AQ 


ec 


XU  Ax 
ec 


AQi  = 


M 


I  2nRTj^  ) 


1/2 


(P.  - 


CpiPiAiATiAx 


At 


^s  = 


Cn  P  A  AT.  Ax 
Ps  s  s  b 

At 


AQ^^  +  AQ,  +  AQ^  =  0  . 
ec  i  s 


By  assuming  no  temperature  drop  at  the  liquid/case  boundary,  i.e.,  Tj^  » 
Eqs.  (19-22)  can  also  be  solved  using  an  iterative  technique  to  obtain  AQ^^ 
AQi,  AQg,  and  AT^^. 

2.3  Condenser  Region 

In  the  condenser  region,  the  boundary  temperature  of  the  heat  pipe  was 


(19) 


(20) 


(21) 

(22) 


assumed  to  be  constant,  resulting  in  governing  equations  similar  to  those  used 
in  the  adiabatic  section. 


AQ  =  XW  Ax 
ec  ec 


I  2nRT 


1/2 


1  ^ 


AQ,  = 


= 


Cp, p,A, AT, Ax 
At 


Cn  p  A  AT.  Ax 
Ps  s  s  b 

At 


AQ  +AQ,  +  AQ=Q^ 
ec  1  s  out 


(23) 


(24) 


(25) 


(26) 


Although  in  normal  operation,  the  condenser  boundary  temperature  would  vary  with 
the  input  power,  in  a  controlled  laboratory  experiment  this  condition  is  easily 
obtained  by  controlling  the  temperature  and/or  flowrate  of  the  condenser 
coolant. 

2.4  General  Assembly 

Because  the  evaporation-condensation  rate  is  extremely  sensitive  to  changes 
in  the  vapor  pressure,  a  predictor-corrector  treatment  was  applied  to  obtain  a 
stable  solution.  Initially,  the  evaporation  rate  was  assumed  to  be  equal  to  the 
condensation  rate  and  the  vapor  pressure  was  calculated.  This  vapor  pressure 
was  then  used  to  compute  the  change  in  the  mass  of  vapor  and  hence,  the 
difference  between  evaporation  and  condensation.  The  original  assumption  of 
equal  evaporation  and  condensation  rates  was  then  checked  and  the  values 
updated.  Preliminary  tests  indicated  that  the  solution  was  stable  regardless  of 
the  time  step  size  used.  In  comparison,  prior  to  using  this  predictor-corrector 
treatment,  the  solution  vibrated  with  increasing  magnitude  for  time  step  sizes 
below  0.01  millisecond. 

The  solution  procedure  consisted  of  initializing  the  variables  and  arrays, 
and  by  iteration,  computing  the  boundary,  liquid,  and  vapor  temperatures  along 
with  the  evaporation  and  condensation  rates,  the  pressure  drops,  and  the 
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flowrates,  for  each  section.  The  resulting  radii  of  curvature  for  the 
liquid-vapor  interface  was  then  computed  and  the  liquid-vapor  mass  flowrate 
corrected  using  the  resulting  vapor  pressure. 

2.5  Transient  Modeling  Results 

Using  the  previously  described  model,  the  cross-sectional  shape  and  size  of 
a  micro  heat  pipe  designed  for  applications  involving  the  thermal  control  of 
ceramic  chip  carriers  were  optimized.  Figure  1  illustrates  this  application  and 
the  general  heat  pipe  configuration.  Figure  2  illustrates  the  critical 
cross-sectional  dimensions  and  shape.  An  interative  approach  was  utilized  to 
develop  the  optimum  shape.  Once  this  had  been  accomplished,  a  free  surface 
energy  minimization  technique^^  was  used  to  determine  the  stable  cross-sectional 
liquid  configuration.  Results  of  this  analysis  indicated  that  regardless  of  the 
orientation  of  the  heat  pipe,  the  majority  of  the  liquid  was  contained  in  the 
upper  corners  (acute  angle)  as  shown  in  Fig.  2,  and  that  only  a  very  small 
fraction  was  contained  in  the  lower  corner  regions. 

Once  the  optimum  shape  had  been  determined,  the  numerical  model  was  used  to 
predict  the  axial  temperature  distribution  of  the  heat  pipe  for  a  given 
operating  condition.  Figure  3  illustrates  the  axial  temperature  distribution 
for  an  input  power  of  0.18  W.  The  temperature  in  the  evaporator  region  shown  on 
the  left  end  decreases  rapidly  with  a  small  continued  decrease  occurring  in  the 
adiabatic  section  and  an  almost  constant  temperature  in  the  condenser  section. 
Three  temperatures  are  indicated;  the  boundary  and  liquid  temperatures  which 
are,  because  of  their  similarity,  almost  identical,  and  the  base  temperature. 
Although  significant  axial  temperature  differences  exist  in  the  boundary  and 
liquid  profiles,  it  is  interesting  to  note  that  vapor  temperature  is  almost 
constant  throughout  both  the  adiabatic  and  condenser  regions. 
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Figure  4  illustrates  the  liquid  mass  flow  rate  as  a  function  of  time  and 
axial  position,  as  was  the  case  for  the  previously  analyzed  trapezoidal  heat 
pipe.  The  tapered  micro  heat  pipe  evaluated  here  displayed  reverse  liquid  flow 
during  startup.  This  reverse  liquid  flow  is  believed  to  be  the  result  of  an 
imbalance  in  the  overall  pressure  drop.  Because  the  evaporation  rate  does  not 
provide  a  sufficient  change  in  the  liquid-vapor  interfacial  curvature  adequate 
to  compensate  for  the  total  pressure  drop  resulting  from  the  temperature  change, 
some  liquid  is  forced  from  the  evaporator  region  towards  the  condenser  region 
through  the  corner  regions  of  the  heat  pipe.  Once  the  heat  pipe  reaches  full 
power,  the  reverse  liquid  flow  disappears,  and  the  liquid  mass  flow  rate 
gradually  increases  until  a  steady-state  condition  is  reached  and  the  liquid 
mass  flow  rate  is  equal  to  the  vapor  mass  flow  rate  in  any  given  section.  This 
phenomena  is  more  apparent  in  Fig.  5,  which  illustrates  the  pressure  as  a 
function  of  axial  position  for  10  second  intervals  approaching  steady-state.  As 
shown,  10  seconds  after  startup  the  pressure  of  both  the  liquid  and  vapor  are 
higher  in  the  evaporator  and  gradually  decrease  with  position.  This  is  the 
phenomena  which  promotes  flow  away  from  the  evaporator  and  towards  the 
condenser.  Once  steady-state  has  been  achieved,  the  liquid  pressure  is  much 
lower  in  the  evaporator  region  and  increases  dramatically  with  position, 
promoting  flow  back  towards  the  evaporator. 


3.0  EXPERIMENTAL  PROGRAM 

To  verify  the  accuracy  of  the  transient  model  and  determine  the  maximum  heat 
transport  capacity,  further  experimental  investigations  were  conducted  using  a 
tapered  micro  heat  pipe  designed  specifically  for  use  in  the  thermal  control  of 
ceramic  chip  carriers.  In  the  application  illustrated  in  Fig.  1,  this  micro 
heat  pipe  fits  securely  under  the  chip  and  is  attached  to  the  chip  carrier.  4 
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Pig.  4  Liquid  nass  flow  rate  as  a  function  of  ti«e  and  aucial  position. 


Small  fins  located  at  the  condenser  end  of  the  heat  pipe  promote  the  cooling  of 
the  condenser  by  free  convection. 

3.1  Experimental  Test  Facility 

The  experimental  test  facility  used  for  this  investigation  has  been 

9 

previously  described  by  Peterson  et  al.  Figure  6  illustrates  the  vacuum  test 

facility  while  Fig.  7  illustrates  the  actual  heat  pipe  test  stand.  Although  the 

test  facility  and  procedure  were  similar  to  the  one  used  for  the  trapezoidal 

heat  pipe,  in  the  present  investigation,  the  heat  source  used  was  a  thin 

nichrome  strip  heater  located  only  on  the  top  portion  of  the  heat  pipe. 

Further,  in  addition  to  a  series  of  eight  copper-cons tantan  thermocouples 

attached  to  the  underside  of  each  test  pipe,  three  in  the  evaporator  section, 

four  in  the  adiabatic  section,  and  one  at  the  inlet  to  the  condenser,  an 

Infrared  thermal  measurement  system  was  used  to  monitor  the  temperature 

distribution  in  the  heat  pipe.  A  Hughes  Probe  Eye  TVS  model  3000  Infrared 

camera  in  conjunction  with  an  RGB  monitor  and  a  real  time  recorder  were  used  to 

observe  the  test  pipe  during  both  startup  and  transient  operation.  Using  this 

system  the  temperature  distribution  along  the  heat  pipe  could  be  measured  with 

an  accuracy  of  +0.05®C.  To  reduce  heat  losses  through  convection,  all  of  the 

_2 

tests  were  conducted  in  a  vacuum  of  <  10  Torr.  The  entire  test  facility  was 
encased  in  a  stainless  steel  vacuum  chamber  equipped  with  a  single  crystal 
silicon  window  which  allowed  observation  with  the  infrared  thermal  measurement 
system.  A  mechanical  roughing  pump  in  series  with  an  oil  diffusion  pump  was 
used  to  obtain  the  required  vacuum. 

3.2  Experimental  Procedure 

In  order  to  fully  understand  the  operational  characteristics  and  performance 
limitations  of  this  tapered  micro  heat  pipe,  an  experimental  test  plan  was 
developed.  This  test  plan  was  comprised  of  three  principal  tasks: 
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Fig.  7  Experimental  Test  Stand 


i.  testing  of  a  charged  silver  heat  pipe  to  determine  the  transient 
characteristics. 

ii.  testing  of  a  charged  silver  heat  pipe  to  determine  the  steady-state 
behavior  for  different  temperatures  and  tilt  angles. 

iii.  testing  of  an  uncharged  silver  heat  pipe  for  comparison  with  the 
charged  heat  pipe  results. 

The  heat  pipe  was  first  inserted  into  the  cooling  chamber  with  the  condenser 
portion  inside.  Then  a  nichrome  resistance  heater  was  attached  to  the  upper 
surface  of  the  evaporator.  Both  the  condenser  and  the  evaporator  portion  were 
1.27  cm  in  length  with  the  remaining  portion  serving  as  the  adiabatic 
region.  The  entire  assembly  was  then  mounted  on  the  test  fixture  as  previously 

4 

described.  Once  the  desired  vacuum  has  been  obtained,  the  cooling  bath 
temperature  was  adjusted  and  the  power  to  the  evaporator  portion  of  the  heat 
pipe  was  turned  on,  thermocouple  measurements  were  taken  at  approximately  0.5 
second  intervals  with  temperature  readings  for  each  of  the  nine  thermocouples 
recorded  as  a  function  of  time.  The  heat  input  for  the  steady-state  tests  was 
increased  in  small  increments,  and  adequate  time  was  allowed  for  each 
incremental  increase  to  allow  the  heat  pipe  to  stabilize.  Once  the  desired  test 
condition  had  been  reached,  temperature  measurements  were  taken  for  a  period  of 
several  minutes  to  insure  steady-state  conditions.  To  obtain  data  for  the  next 
successive  power  level,  the  temperature  of  the  cooling  bath  was  adjusted  to  hold 
the  temperature  at  the  entrance  to  the  condenser  region  constant. 

3.3  Experimental  Results 

Figure  8  illustrates  transient  temperature  characteristics  of  the  tapered 
micro  heat  pipe  for  an  input  power  of  0.508  W  and  a  condenser  temperature  of 
23“C.  The  abscissa  is  shown  as  a  temperature  difference  between  the 
thermocouple  located  at  the  inlet  to  the  condenser  and  the  one  indicated  in  the 
legend.  As  illustrated,  the  heat  pipe  reaches  a  steady-state  value  very  quickly 
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Pig.  8  Transient  temperature  distribution  (Power  0.058  V,  T  23”C). 


with  approximately  90%  of  the  maximum  power  obtained  by  approximately  20 
seconds.  At  a  time  less  than  20  seconds,  Thermocouples  1  and  3  yield 
approximately  identical  readings.  At  15  seconds,  however,  they  begin  to 
separate  and  the  temperatures  start  to  diverge.  This  rapid  separation  is 
believed  to  be  the  result  of  the  onset  of  dryout.  Note  that  both  Thermocouples 
1  and  3  are  considerably  hotter  than  Thermocouple  2.  This  may  be  due  to  the 
presence  of  the  lead  wires  at  each  end  of  the  heater  but  this  justification  was 
not  verified.  When  there  is  no  dryout  of  the  corner  regions,  the  temperature 
drops  between  Thermocouples  1  and  2,  and  2  and  3,  respectively,  will  be 
approximately  equal. 

Figure  9  illustrates  the  temperature  difference  between  each  of  the 
thermocouples  as  a  function  of  the  power  input.  Again  the  abscissa  is  shown  as 
a  temperature  difference  between  the  thermocouple  located  at  the  inlet  to  the 
condenser  and  the  one  indicated  in  the  legend.  As  illustrated,  the  heat  pipe 
remains  approximately  isothermal  along  the  entire  length  for  powers  between 
approximately  0.5  and  1.3  W.  Prior  to  the  power  level  of  0.5  W,  a  majority  of 
the  heat  is  either  conducted  axially  through  the  wall  of  the  heat  pipe  or 
radiated  to  the  environment.  At  a  value  of  1.3  W,  rapid  increases  in  the 
evaporator  temperature  (as  Indicated  by  Thermocouples  1,  2,  and  3)  occur.  These 
increases  cascade  throughout  the  evaporator  indicating  the  dryout  has  begun  and 
is  progressing  along  the  length.  This  is  more  apparent  in  Fig.  10,  which 
illustrates  the  same  experimental  results;  however,  the  axis  has  been  enlarged 
to  clarify  the  temperature  differences.  At  approximately  2  W,  this  dryout  has 
reached  Thermocouple  4  located  in  the  adiabatic  section.  This  dryout  "front" 
then  proceeds  reaching  Thermocouple  5  at  a  power  of  approximately  3  W.  It  is 
interesting  to  note  that  at  approximately  0.4  W,  Thermocouples  3  and  4  invert  as 
do  Thermocouples  6  and  7  at  a  power  of  approximately  0.9  W.  No  explanation  has 
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Fig.  10  Enlarged  steady-state  temperature  distribution  as  a  function  of  the 
input  power  (No  tilt,  T  24.5*C). 


been  found  to  determine  the  cause  of  this  phenomena,  however,  it  was  shown  to  be 
repeatable. 

As  mentioned  previously,  at  power  levels  below  0.5  W,  much  of  the  heat  is 
given  off  by  radiation  or  conducted  axially  through  the  wall  of  the  heat  pipe. 
This  is  more  clearly  apparent  in  Fig.  11,  which  illustrates  the  experimental 
*  results  for  an  uncharged  heat  pipe  using  the  same  condenser  temperature  and 

power  levels  used  for  the  previous  figures.  As  illustrated,  the  temperature 
difference  indicated  by  the  straight  line,  has  the  same  slope  for  each  of  the 
corresponding  thermocouples  indicating  that  at  these  power  levels,  the  axial 
temperature  drop  is  the  same  for  both  the  uncharged  and  charged  heat  pipes. 
However,  as  the  power  level  increases,  the  heat  pipe  begins  to  operate  more 
efficiently,  and  the  axial  temperature  drop  remains  constant.  This  phenomena  is 
similar  to  that  observed  in  very  short  heat  pipes.  At  low  powers  they  are 
relatively  inefficient  owing  to  the  thermal  resistance  of  the  case,  wick,  etc., 
but  as  the  length  increases,  the  advantages  become  clear. 

To  determine  the  effect  of  tilt  angle  on  the  operation  of  the  tapered  micro 
heat  pipe,  several  tests  were  conducted  with  similar  input  powers,  but  at 
different  tilt  angles.  Figure  12  illustrates  the  experimental  results  for  the 
tapered  heat  pipe  horizontal  orientation  at  a  power  of  0.12  W  and  a  condenser 
temperature  of  25‘’C.  Figure  13  shows  an  orientation  15  degrees  with  the 
•i  horizontal  (evaporator  up),  and  Figure  14,  at  an  orientation  of  15  degrees  from 

the  horizontal  with  the  evaporator  down.  As  illustrated,  the  experimental 
results  demonstrate  very  little  difference  in  the  temperature  distribution  along 
the  pipe  as  a  function  of  the  input  power  because  of  tilt. 

3.4  Comparison  of  Experimental  and  Analytical  Results 

Figure  15  compares  the  experimental  data  with  the  temperature  distribution 
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Fig.  12  Transient  behavior  prior  to  the  onset  of  dryout,  (Power  0.12  W, 
24. 5*0,  no  tilt). 
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Pig.  13  Transient  behavior  prior  to  the  onset  of  dryout,  (Power  0.12  W, 
24.5"C,  15  degree  tilt  evaporator  up). 
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Fig.  14  Transient  behavior  prior  to  the  onset  of  dryout,  (Power  0.12 
24.5®C,  15  degree  tilt  evaporator  down). 
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Fig.  15  Coaparison  of  predicted  and  measured  boundary  temperature  distribution, 
(Power  0.12  W,  T  24.5"C,  no  tilt). 


as  predicted  by  the  numerical  model  for  a  power  level  of  0.12  W.  As 
illustrated,  the  temperature  data  closely  follow  the  predicted  trend.  Some 
heat  rejection  occurrs  in  the  adiabatic  section,  however,  the  condenser 
temperature  remains  relatively  constant.  Throughout  the  entire  length,  however, 
the  predicted  and  measured  temperature  distributions  are  extremely  close  to  the 
r  predicted  value,  with  the  largest  deviation  being  less  than  0.3'’C.  Figure  16 

compares  the  predicted  and  the  measured  temperature  difference  between  the 
evaporator.  Thermocouple  1,  and  the  condenser.  Thermocouple  9,  as  a  function  of 
time  for  a  power  input  of  0.12  W  and  at  a  horizontal  position.  As  illustrated, 
the  final  steady-state  temperature  difference  between  the  evaporator  and  the 
condenser  as  predicted  by  the  model  is  within  12Z  of  the  measured  value.  The 
numerical  model,  however,  predicts  a  much  more  rapid  startup  of  the  heat  pipe 
than  that  obtained  experimentally.  This  predicted  rapid  startup  may  be  the 
result  of  the  Inherent  shortcomings  of  the  model,  while  in  the  actual  case  heat 
is  rejected  from  the  evaporator  portion  of  the  heat  pipe  through  vaporization  of 
the  liquid,  it  is  also  rejected  through  conduction  to  the  support  member  and 
radiation  to  the  surroundings.  In  addition,  the  numerical  model  assumes 
instantaneous  heat  addition  while  in  the  experimental  case,  some  sensible  heat 
is  absorbed  by  the  mass  of  the  heater  wires,  the  heater,  and  the  lead  wires. 

The  small  dip  in  the  theoretical  predicted  temperature  difference  occurring  at 
f  approximately  10  seconds  into  the  test  is  the  result  of  the  reverse  liquid  flow 
discussed  previously.  As  shown,  the  experimental  results  indicate  no  such 

t 

reverse  flow.  This  may  be  the  result  of  the  reduced  time  constant  caused  by  the 
mass  of  heat  pipe  case  material. 

Figure  17  illustrates  a  similar  comparision  at  a  tilt  angle  of  15  degrees 
(evaporator  down).  As  shown,  for  this  case,  the  numerical  model  does  not 
predict  any  reverse  liquid  flow  and  significantly  underestimates  the 
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Fig.  16  Conparison  of  predicted  and  measured  evaporator  temperature  difference 
as  a  function  of  time,  (Power  0.12  W,  T  24.5“C,  evaporator  down). 
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^^8"  Conparlson  of  predicted  and  measured  evaporator  temperature  difference 
as  a  function  of  time,  (Power  0.12  W,  T  24.5"C,  no  tilt). 
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longitudinal  temperature  difference.  Again,  the  disparity  between  the  measured 
and  predicted  values  is  believed  to  be  due  to  differences  in  the  numerical  model 
and  the  actual  experimental  tests.  In  the  numerical  model,  the  effects  of  the 
end  wall  of  the  heat  pipe  were  not  incorporated.  This  end  region  provides  an 
additional  capilllary  radius  which  helps  retain  liquid  in  the  evaporator  region 
of  the  heat  pipe  and  results  in  higher  power  levels  prior  to  the  onset  of  ' 

dryout.  These  increased  power  levels  manifest  themselves  as  an  increase  in 
temperature  differences,  similar  to  that  shown  in  Fig.  17. 

4.0  SUMMARY 

An  investigation  to  develop  a  numerical  model  capable  of  predicting  the 
transient  performance,  determining  the  effect  of  variations  in  the  size  and 
optimizing  the  cross-sectional  shape  of  a  tapered  micro  heat  pipe  has  been 
conducted.  Several  heat  pipes  which  utilize  the  optimized  shape  were  evaluated 
experimentally  to  determine  their  suitability  for  heat  removal  from  a  ceramic 
chip  carrier  and  to  provide  data  through  which  the  numerical  model  could  be 
validated.  The  numerical  model  was  found  to  be  capable  of  accurately  predicting 
the  maximum  transport  capacity  prior  to  the  onset  of  dryout,  the  temperature 
distribution  throughout  the  longitudinal  position,  and  the  temperature 
difference  between  different  locations  on  the  heat  pipe  to  within  prior  to 

the  onset  of  dryout.  In  addition,  the  transient  model  predicted  the  operation  f 

and  temperature  distribution  of  the  tapered  micro  heat  pipe  to  within  ±12  X  over 
a  wide  range  of  powers  and  operating  temperatures. 

As  a  result  of  this  effort,  the  numerical  model  can  be  confidently  used  to 
identify,  evaluate,  and  better  understand  the  phenomena  that  governs  the 
transient  behavior  of  micro  heat  pipes  as  a  function  of  the  physical  shape, 
properties  of  the  working  fluid,  and  principal  dimensions. 
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